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Abstract 

Quarkonium production in high-energy reactions is found to exhibit a be- 
haviour more universal than that expected from velocity scaling. Total 
rates of quarkonia produced in hadronic interactions as well as Feynman- 
x and transverse momentum distributions can be described over the full 
range of accessible energies (15 1800 GeV) by two-stage processes. 

The quarkonium production cross section factors into a process-dependent 
short-distance part and a single long-distance matrix element. The first 
part describing the production of a free quark-antiquark pair is the pertur- 
batively calculated subthreshold cross section. The non-perturbative factor 
turns out to be universal, giving the model great predictive power. Further- 
more we estimate the fraction of the heavy-quark cross section leading to 
quarkonium for both the charm and bottom systems. Finally, we comment 
on quarkonium photoproduction. 
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Recently much effort has been devoted to explain quarkonium production in a new 
factorization approach [|l| . Any quarkonium cross section is given by the sum of infinitely 
many terms, each of which factors into the product of two terms. The first one, calculable 
as a series in a s (n) where /i is of the order of the relevant hard scale, is the cross section 
to produce a free quark-antiquark pair (QQ) in a particular angular momentum and 
colour state 2S + 1 l,y (in the spectroscopic notation and c = 1 (c = 8) denotes a colour- 
singlet (colour-octet) state). The second factor determines the probability that such a QQ 
pair binds to form a quarkonium H{n.J PC ) of given total spin J, parity P, and charge 
conjugation C . The factorization approach becomes meaningful through the velocity- 
scaling rules, which determine the relative importance of the various long-distance matrix 
elements (ME). At any order in v, the velocity of the heavy quark within the bound state, 
the quarkonium cross section is hence given by a finite number of contributions @. 

This results in an expansion of the quarkonium cross section in both a s (fj,) and v. Al- 
though this "velocity-scaling model" (VSM) suggests an explanation of quarkonium pro- 
duction at the Tevatron , its weak point is the fact that currently the non-perturbative 
ME cannot be calculated in QCD[|. A crucial test of the approach is therefore the de- 
termination of the MEs from as many different high-energy reactions as possible. This 
endeavour is, however, rendered more difficult by the fact that, in general, different combi- 
nations of MEs arise. Nonetheless, preliminary attempts indicate that the velocity-scaling 
is not perfect: J/ip production at the Tevatron requires considerably larger c = 8 MEs 0] 
than J/ip hadropro duct ion || or the ^-distribution in J/ip photoproduction ||. Also the 
hadroproduction ratio Xci/Xc2 is too low compared to data since Xci production is clearly 
disfavoured by either a power of a s (m c ) or a factor (v 2 ) 2 compared to Xc2- Last but not 
least, the J/ip (non-) polarization is difficult to account for in the VSM |5], [7J. 

Some time ago we || have shown that existing quarkonium production data at fixed- 
target energies are, in fact, compatible with the assumption that the non-perturbative 
transition of the QQ pair to quarkonium is more universal than expected from the velocity- 
scaling rules. Indeed, low-energy data are well reproduced if the (infinite) sum of short- 
distance coefficients times long-distance MEs is truncated to a single term 



a 



H (nJ PC 



F[nJ PC } & QQ . (1) 



In this letter we shall demonstrate that the colour-evaporation model (CEM) of eq. ([!]) also 
accounts for quarkonium production at the Tevatron and comment upon its application 
to photoproducton of quarkonia. We emphasize that in contrast to the VSM, the number 
of non-perturbative parameters is minimal, hence the model possesses great predictive 
power. 

The short-distance part of eq. (P is the perturbative subthreshold cross section ex- 
panded in powers of a s (/i) where \i oc uiq. Specifying to charm, the cross section is 

where cx[cc] is the spin- and colour-averaged open heavy-quark pair production cross sec- 
tion. The CEM is hence based on two ingredients. First, the quarkonium dynamics are 
assumed to be identical to those of low mass open QQ pairs. All perturbative QCD 
corrections apply to the short-distance cross section and hence are identical for open 



1 Attempts to calculate decay MEs on the lattice have just started 
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Figure 1: Left: The J ftp production cross section cr pN [J/?/>] for xp > 0, calculated 
with MRS D-' PDF, compared to data |J. Right: Energy dependence of T production 
y u+ / u"]da[T(nS)]P N /dy at y 



X] n Br[T(nS) — > ]da[T(nS)] piN /dy at y — compared to data [T(], [nj; the predic- 
tions with MRS D-' and GRV HO PDF essentially coincide. Also shown (CR) is the 



phenomenological fit of [12|; from 
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Figure 2: The J/ip longitudinal momentum distributions compared to piV and pN data 
l|, with x F = p L [J/^]/Pmax[J/ip]; results obtained with the MRS D-' (GRV) PDF are 



denoted by a solid (dashed) line; from 
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and bound heavy-quark production. Second, although the QQ pair is produced at short 
distances in different states (distinguished by colour, angular momentum, relative mo- 
mentum) and their relative production rates may (and will) be different for different 
high-energy collisions, it is only the average over many long-distance matrix elements, 
combined in the long-distance factor F[nJ }, that determines the probability to form 
a specific bound state. Necessarily, the factor F needs to be universal, i.e. process- and 
kinematics- independent. 

To illustrate the success of the CEM, in Fig. [I] we compare the prediction for the total 
J I ip and T production rates with data. Note that the model uniquely predicts the shape 
of the energy dependence while the absolute normalization at low energies fixes the non- 
perturbative factor F. Fig. || shows the prediction of fixed-target xp distributions. There 
is remarkable agreement over a wide energy range, from low-energy pp collisions where 
valence qq fusion dominates up to high-energy pp collisions dominated by gluon-gluon 
fusion. Note that both the shape and normalization of the xp distributions are obtained 
from the model once F has been fixed by the energy dependence. 

The long-distance factors determined from the low-energy total cross sections in Figs. p] 

are 



F tat [J/i/>] = 2.5% 
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J2 Br[T(nS) -> /i VI F tot [T(nS)\ = 1.6 x 1(T 3 . (3) 



n=l 



Here the subthreshold cross sections were calculated in next-to-leading order (NLO) using 



the MRS D-' parametrization JM] of the parton-distribution functions (PDF) with renor- 
malization and factorization scales fiR and fip chosen to be [Ir = fip = 2m c = 2.4 GeV 
and fiR = fip = m b = 4.75 GeV, respectively^. The results in eq. (f|) are rather in- 
sensitive to variations of the parameters in the open heavy-quark cross section, if they 
are tuned to the open heavy-flavour total cross section data. For instance, the GRV 



HO parametrization JT5J with Hr = Hf = m c — 1-3 GeV leads to very similar results: 
the smaller subthreshold region is basically compensated by the larger two- loop a s value, 
0.298 for GRV HO (A 4 = 0.2 GeV) compared to 0.243 for MRS D-' (A 4 = 0.23 GeV). Note, 
however, that the long-distance factors will be considerably larger if the open heavy-quark 
cross section is calculated to leading order (LO) only, a factor 2.2 larger for J/ip and a 
factor 1.9 larger for T. 

In Fig. |H we show the transverse momentum distributions of prompt charmonium 
production (i.e. not originating from B decays) at the Tevatron energy y^i =1.8 TeV. The 
bottomonium transverse momentum distributions are given in Fig. |j. The normalizations 
for the various states are given in table [I], obtained from eq. using the measured cross 
section ratios || and the branching ratio to muon pairs fL8"f . Good agreement with data 
is found, typically better than 50%. Note that the CEM prediction shown in figs. |3] and |] 
is based on the subthreshold cross section calculated to O(a^), which is NLO for the pp 
integrated cross section but LO only for the pp distributionP]. In the absence of the NLO 
corrections to the pp spectrum one might apply a K factor to account for the unknown 
higher-order corrections or simply use F-values extracted in LO (two such examples are 

2 To calculate the NLO subthreshold cross section, we use the program of Mangano, Nason, and Ridolfi 
p6[ , restricting the QQ mass range. 

3 In the calculation of the pp distribution, we used fjr R = n 2 F — n 2 [niQ + {pp q +p^^)/2] with n = 2 

for charm and 1 for bottom consistent with the scale advocated for open QQ production |16|. 
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Figure 3: Transverse momentum distributions of charmonia compared to CDF data [T7[ . 
The upper curve for ip(2S) contains an extra i^-factor of 2. 
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Figure 4: Transverse momentum distributions of bottomonia compared to CDF data 
| n| for various values of the intrinsic transverse momentum with F as in table 1: (hx) = 
1.25 GeV (dotted), 2.0 GeV (solid), 2.5 GeV (dot-dashed). The T(1S) prediction without 
smearing {kj) = GeV is shown as dashed lines, the upper curve containing an extra 
K-i actor of 1.9. g 
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Table 1: Long-distance factors F tot from eq. (D) expressed in percent for 1 states and 
the sum of the inclusive Xcj production rates (i.e. including cascade decays). The factors 
which multiply the NLO subthreshold cross sections do not include the branching ratio 
into lepton pairs. For LO calculations the above numbers should be multiplied by 2.2 for 
charmonia and 1.9 for bottomonia. 



given in Figs. ||] and §). In either case one expects larger cross sections so that our 
estimates are rather conservative. 

The CEM prediction for ip(2S) is about a factor of two lower than the data. This 
is a simple consequence of the fact that at high the ip(2S)-to-J/ip ratio measured 
at the Tevatron is about twice as large as that observed at fixed-target energies for the 
Pr-integrated cross section. Multiplying the CEM prediction by a factor of two produces 
very good agreement. More precise data will show whether this is a systematic effect that 
would require refinements of the CEM. 

The CEM prediction for the pt distributions based on fixed-order perturbation the- 
ory cannot be trusted for ^rriQ. A correct treatment of the \ow-pr region requires 
soft-gluon resummation and the inclusion of intrinsic transverse momenta, analogous to 
the Drell-Yan case. The effect of soft-gluon resummation can be mimicked through an 
effective, larger value of the average intrinsic transverse momentum (kj-). Fig. |] shows 
that inclusion of (kx) smearing results in good agreement with data down to very low pr- 

The information of table can be used to estimate the total bound-state probability. 
In the case of charmonium this requires assumptions about the rj c (nS) and h c (lP) cross 
sections. In the case of bottomonium, we also need assumptions about the ratios of the 
Xbj(^-P) to T(nS) cross sections. Taking the latter to be equal to that measured in the 
charmonium system and assuming that both the S- and P-wave state cross sections are 
proportional to 2 J + 1 and disregarding possible bb .D-wave states we estimate 

F[i] ~ (8 - 10)% charm 

i 

« (17 -32)% bottom . (4) 

The dominant part of the subthreshold charm cross section produces open charm. This 
fraction is considerably reduced in the bottom system: We observe a significant increase of 
the bound state fraction with increasing quark mass. The fact that the total charmonium 
cross section is just 1/(1+8) of the subthreshold cross section must therefore be considered 
as fortuitous. This ratio was recently advocated as universal for colour- singlet production 
TT| giving the fraction of both diffractive events in deep- inelastic ep scatterings and bound 



states in heavy-quark production. Our analysis shows that bound state production does 
not obey this rule. Moreover there is no reason to expect this ratio to hold. In fact, 
considering the complete system rather than restricting to the QQ pair suggests that the 
colour-singlet fraction is 1 : 1 rather than 1 : 9 p0 |. 



Finally we discuss photoproduction of charmonium. We first note that the inelastic 
J/ij} cross section is defined only within cuts. These cuts are necessary since the quasi- 
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elastic process, 7p — ► J ftp + p, and the forward-elastic reaction, 7p — > J/-0 + X, where 
the is isolated in rapidity, cannot be described in perturbative QCD using a (single) 
gluon distribution in the nucleon. In order to stay away from diffractive production one 
typically restricts z (E^/ E 1 in the nucleon rest frame) to values less than 0.9 and applies 
additional cuts as, e.g., on p T or on the number of tracks. Since the cross section rises 
quickly towards 2 = 1, F[J ftp] is not well defined heref]. Comparing the cross sections for 
open charm and (inelastic) J / xp photoproduction we find 

F[Jftp) « (1 - 2.5)% at low y/a 

« (0.5 - 1.4)% at high y/s , (5) 

numbers in reasonable agreement with the value obtained from hadropro duct ion. While 
the J / ip long-distance factor thus appears to be universal, photoproduction of Xcj seems 
to be problematic for the CEM: NAM |2T| puts an upper limit of about 8% for the 
fraction of J ftp's from Xcj decays while hadroproduction experiments suggest a value 
four times greater. So far these measurements are not very precise, but, if confirmed, 
Xcj photoproduction might indicate limitations on the CEM, requiring refinements of its 
simplest variant. 

In summary, we find an impressive agreement between data on quarkonium produc- 
tion in hadronic collisions and the CEM: Total cross sections, xp distributions, and pp 
spectra are all well described by the assumption that the non-perturbative bound-state 
formation is governed by an average, universal, long-distance factor. Since only a single 
non-perturbative ingredient is required for any given bound state, the CEM has great pre- 
dictive power. Two possible deviations from this simple scenario have been pointed out: 
a difference in the tp(2S)f J ftp ratio in hadronic collisions and a different Xcj production 
fraction in hadro- and photoproduction. More data, also for bottomonium, are eagerly 
awaited. 

G.S. thanks the SLAC theory group for hospitality during the completion of this work. 
We also thank H. Satz for many fruitful discussions. 
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